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TECHNICAL NOTE 3570 

AN EXPERIMENTAL COMPARISON OP THE LAGRANGIAN MD EUEERIAN CORRELATION 

4 

COEEFICIENTS IN HOMOGENEOUS ISOTROPIC TURBULENCE 
By William R. Mickelsen 


SUMMARY 

An experimental comparison was made “between the Lagrangian and 
Eulerian correlation coefficients in the hcmogeneous, isotropic central 
core of a turbulent pipe flow. The Lagrangian correlation coefficient 
was characterized by measurements of the turbulent diffusion of helium 
from a point source. The turbulence intensity and Eulerian correlation 
coefficients were measured by hot-wire anemometry. 

The Lagrangian and Eulerian correlation coefficients had similar 
shapes related by a linear proportionality between their time and space 
coordinates. The linearity and the proportionality factor were evaluated 
over a range of turbulence intensities from 1.8 to 14 feet per second. 

The proportionality factor was roughly 0.6 and provides a means for 
the quantitative solution of practical turbulent-mixing problems from 
hot-wire -anemometry data or from knowledge of the Eulerian parameters 
of the turbulence field. 


INTRODUCTION 

The tiirbulent mixing of mass and heat is of major importance to the 
performance of jet-engine combustors. For example, mixing is largely 
responsible for the preparation of desired fuel-air mixtures and also 
affects the over-all combustion process. The general objective of this 
investigation >jas to find relations between turbulent -milxing theory and 
experiment which would provide sufficient Information for the direct so- 
lution of practical mixing problems. 

The Taylor theory of diffusion by continuous movements (ref. l) 
provides definitions of the turbulent parameters which enter directly 
into the solution of tirrbulent-mixing problems. This theory will be 
disciissed in detail in the next section of this report, and it will 
sirfflce to say here that the turbulent-mixing parameters eire derived 
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directly frcan tiie LagrangiaQ correlation coefficient 9 . . This coeffi- 
cient expresses the degree of correlation between the velocity v^ of 

a particle at time tQ and the velocity v^ of the same particle at a 
later time t, as shown hy the equation 

V V 

« - ( 1 ) 

(See appendix A for definitions of symbols.) Equation (l) is written 
for a homogeneous field of turbulence, and the bar over the velocity 
products denotes the average of a large number of particle motions. The 
Lagrangian correlation coefficient Si may be further clarified by the 
foUoifing figures; 



Position of particle Time, t, sec 

at zero time 

As shcnm in the right-hand graph, ^ = 1 at t = 0, and 31 usually ap- 
proaches zero as t approaches infinity. 

Although the Lagrangian coiTrelation coefficient cannot be measured 
directly by existing instrumentation, the mixing rate of any particular 
turbulent field can be determined by measuring the diffusion of mass or 
heat from a simple source. The Eulerian turbulence parameters of the 
field can be measured by hot-id.re anemcmetiy, or possibly predicted from 
experience or from the shape of the flow- field boundaries. The direct 
measurement of the mixing rate is laborious and awkward compared with 
hot-'td.re-anemometry techniques j hence, a need exists for knowledge of re- 
lations betireen the Lagrangian and Eulerian parameters of turbulence. 

Two important Eulerian measures of a tiorbulence field are the turbu- 
lence intensity a/? and the double velocity correlation coefficient, a 
nine -component second-order tensor, which for homogeneous isotropic tur- 
bulence can be reduced to two interrelated correlation coefficients f 
and g, as shorn by Kaim^n and Howarth (ref. 2). The Eulerian correlation 
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coefficients f and g express the degree of correlation "between 
tiirhulent-velocity components of different fluid particles at points 
separated by a distance 


u 


o 





u^ 

U Uj. 


€ 

^ O t 

(2) 

1 

u^ 


f 

V Vf 

(3) 

^ — H 

v^ 



where the bars denote time averages. The relation between f and g 
is shown in the next section of the report. Both f and g become 
equal to 1.0 at ^ = Oj and both usually are asymptotic to zero as ^ 
approaches infinity, in a manner similar to that of the Lagrangian cor- 
relation coefficient ^ . Both f and g can be conveniently measured 
by present-day hot-wire anemometry. 


The only work appearing in the literature on the theoretical rela- 
tion between the Lagrangian and Ihilerian correlation coefficients is 
that of reference 3. The analysis is based on Joint-velocity probability 
distributions and Eulerian correlation coefficients between successive 
Incremental distances but is not yet in a completely usable form. 

A mmber of experimental measurements of the mixing of mass and 
heat have been made in turbulent fluid streams. In the work of refer- 
ence 4, hydrogen and carbon dioxide were injected from a simulated point 
source into the central core of a fully developed turbulent, pipe air 
flow and vrere found to diffuse at the same rate. The results agreed with 
Taylor's theory of diffusion in that the turbulent diffusion coefficient 
was a function of time which Increased from zero to an asymptotic value 
far. downstream. 


A similar experiment was performed by measuring the turbulent dlf- 
fixslon of a salt solution from a simulated point source in a fully 
developed turbulent water channel flow (ref. 5). The results, which 
agreed with those of reference 4, showed the Lagrangian correlation co- 
efficient to be a continuoiis function of time having the characteristics 
discussed in a preceding paragraph. 

Measirrements have also been made of the turbulent diffusion of heat 
close behind a line source in a grid-generated flowing turbulence field 
(refs. 6 to 8) . The Taylor theory of diffusion predicts that, for very 
short times subseq.uent to the start of the diffusion process, the turbu- 
lent diffusion coefficient is equal to v^, where v^ is the mean 
square of the turbulent- velocity f luct\aations , and t is the residence 
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■time in the stream. Reference 6 reports work -vAiere the turbulent- 
diffusion coefficient was equal to that predicted, by the Taylor theory, 
while references 7 and 8 give experimental turbulent-diffiision coeffi- 
cients hl^er than those predicted theoretically. 

The experimental work of references 7 and 9 included determination 
of the Lagrangian correlation coefficient in fields of decaying iso- 
tropic turbulence. These measurements are also in accord with the 
Taylor description of the Lagrangian correlation coefficient. 

The specific objective of the present investigation was to find a 
relation bet\reen the Lagrangian and the Eulerian turbulence parameters, 
so that q\iantitative solutions to mixing problems could be obtained 
from anemometry measurements of the turbialent field. Lagrangian 
turbulent-mixing parameters were found experimentally by measuring the 
diffusion of helium from a point source in an air stream. Helium was 
used as a diffusing material because it could be detected in very low 
concentrations with existing mass-spectrometry equipment. The low con- 
centrations were desirable since the helium injection flow rates, and 
hence the helium injector, could be kept small. In this way, the ef- 
fect of the injector vortex wake on the diffusion process was minimized. 

In addition, the loi-r concentrations of helium had a negligible effect 
on the fluid properties of the stream. The Eulerian correlation coeffi- 
cients were measiored -id-th constant-temperature-anemometry equipment. 

The experiment was performed in the approximately homogeneoxis, isotropic 
central core of fully developed turbulent pipe flow. 

The investigation was carried out at the MCA Leid.s Plight Propulsion 
Labonratory as a part of the conibustion research program. 


THEORETICAL ANALYSIS 

Defining the turb-ulent mixing downstream from a point source in terms 
of the Lagrangian correlation coefficient provides a comparison between 
the Lagrangian and Eulerian correlation coefficients. This section shows 
the derivation of this Lagrangian-mixlng relation and also includes 
Eulerian relations pertinent to the condition of isotropy and to the 
Lagrangian-Eulerian conparlson. 


The Taylor Theory of Diffusion 

Taylor shows in reference 1 that the mean displacement of fluid 
particles in a turb-ulent field is related to -the Lagrangian correlation 
coefficient 


1 ^ 
2 dt 



9t dt 


( 4 ) 
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o 

where Y is the mean square displacement of a group of many particles. 
Reference 10 points out that hy analogy to Brownian movement 


1 dY^ 
“ " 2 dt 


(5) 


where D Is a turhulent-dlf fusion coefficient. 

As noted In reference 11^ experiment (refs. 4 to 9) has Indicated 
that In homogeneous turbulence^ the fluid particle displacements Y 
have a Gaussian probability density. Reference 11 shows that as a con- 
sequence of this Gaussian form the diffusion coefficient defined by equa- 
tion (5) may be tised in the Fourier-Poisson diffusion eguation 

If = (i = 1,2,3) (6) 

where C Is the concentration of mass at the coordinates at time 

t, and d is the molecular-diffusion coefficient. 


Solution for the Point Source of Diffusion in a Moving Stream 


If the molecular-diffizsion coefficient d is constant, as in the 
present investigation, and if the tvirbxilent field is homogeneous and 
isotropic, then equation (6) can be written 


If = = 1^2,3) 

A new parameter O) may be defined by the relation 


0) = oij + CD^ = (D + d) dt 

vJo 


(7) 


( 8 ) 


where £ 0 , oOp, and are called the total, the turbulent, and the 

molecular-spreading coefficients, respectively. Through this transforma- 
tion, equation (7) may be written 


3C ^ ^ 
. Sx| 


(1 = 1,2,3) 


(9) 
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The solution to equation (9) for an instantaineous point source re- 
leased into a stream having a velocity U is given in reference 12 as 


C = 


W 


(Aitm) 




(x-Ut)^+r^ 

4o) 


( 10 ) 


Reference 13 notes that the solution for a continuous point source can 
he obtained hy summing the contributions from the infinite number of 
diffusion centers downstream from the injection point; 




c = 


w 


(4«) 


3y^ 




(D 


3/2 


_ (x-Ut) ^+r^ 
40) 


dt 


( 11 ) 


An approximate solution for the continuous point source suggested 
in reference 4 for cases where the diffusion in the x direction is 
negligible and where x » r is 


C =_5L_ 
43tttD 



( 12 ) 


The validity of equation (12) was checked in the present investigation 
by comparing it numerically -^d.th a graphical solution of equation (U) . 
The comparison indicated that equations (u) and (l2) are practical .ly 
identical over a range of r sufficient for the evaluation of cd. 


Relations Between the Eulerian Correlation Coefficients in 
Homogeneous Isotropic Turbulence 

For homogeneoTis, Isotropic turbulence, Karm^ and Howarth (ref. 2) 
show that the Eulerian double velocity correlation coefficient is a 
second-order tensor, that is. 


f 0 0 
0 f 0 
0 0 g 


(13) 


where the correlation coefficients f and g (eqs. (2) and (3)) are re- 
lated through the equation 

2 d^ 


f - g 


(14) 
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Another requirement for the condition of Isotropy is that 

(l5) 

Equations (l4) and (l5) provide means for ascertaining the degree of 
isotropy in the experimental turbulent stream. 


Comparison of the Lagrangian and Eulerian Correlation Coefficients 
A Lagrangian scale of turbulence jSf may be defined by the relation 



t* 


where t* is a characteristic time given by 

X OO 

^ dt 


Combining equations (l6) and (l7) gives 


if = 



dt 


In an analogous manner, an Eulerian scale can be defined by 

J ’^0® 

0 


(16) 


(17) 


(18) 


(19) 


For a homogeneous field of turbulence, the Lagrangian and Eulerian 
scales are proportlonalj 


Lf = Bif (20) 

Furthermore, if the Lagrangian and Eulerian correlation coefficients 
have the same functional forms, then a basis for their ccaiparison can 
be expressed by 


f(^) 



( 21 ) 


where B is a numerical factor of proportionality. Introduction of 
equations (20) and (21) into (l9) yields equation (18), as consistency 
requires . 
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The Lagrangiaii data taken in the present investigation consisted of 
the spreading coefficient cd whi.ch, when corrected for the molecular 
spreading throu^ equation (8), was reduced to the t\irhulent-spreading 
coefficient oip. Ccmhi nation of equations (4)^ (5), and (8) results in 

the following escpression for 

— pt pt 

= v^ J J « dt dt (22) 

Rather than perform a douhle differentiation of the data, it was 

deemed more accurate to perform a douhle integration of the Eulerian 
correlation coefficient f hy the relation 

By comparison between the cuf calculated from the anemometry data and 
the obtained from the helivtm-diffusion data, the similarity hypothe- 

sis expressed by equation (2l) could be Investigated. 

As shotm by the discvission in this section, it is possible to com- 
pare the shape characteristics of the Lagrangian and Eulerian correlation 
coefficients by a combination of helium-diffusion and hot-t'dre-anemometry 
measurements . 


APPARATUS 

Ducting and Air-Flow Instrumentation 

The turbulent field ^vas generated in an 8-inch-diameter duct having 
an inlet length- to-diameter ratio of .36. As shown in figure 1, room air 
was dra-vm into the duct through an inlet contraction nozzle and twO' 
screens to promote nearly fully developed pipe turbulence at the helium 
injector and instrumentation stations. The air flow passed throu^ a 
critical-flow "exiaTist choke" station and then into the laboratory alti- 
tude exhaust facility. The air-flow rate was controlled by varying the 
exhaust eirea. 

The air temperature was measvired with a mercury thermometer located 
outside the duct near the inlet nozzle. Duct static pressures vrere meas- 
ured -idth standard 1 taps placed at the Instrumentation stations and i 

connected to mercury manometers. 

The stream velocity "tTas measured with a standard pitot static probe 
connected to a micromanometer and with a hot-wire anemometer by using a 
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vortex- shedding technique described in reference 14. The techn i que con- 
sists of determining the frequency of vortex- shedding from a cylinder 
with a hot-v/ire anemometer. The frequency of vortex-shedding n is re- 
lated to the stream velocity U by the StrotJhal number S' = ndg/lJ, 

where the Strouhal number S is constant over a large range of Reynolds 
number based on the cylinder diameter d^. The hot-wire-anemometry 
equipment also provided a convenient means for measuring radial and 
axial stream- velocity profiles. 


Turbulence Instrumentation 

Hot-wire-anemometry equipment was used to measure the turbulence 

intensities and and the Eulerian correlation coefficients f 

and g. The hot-wire anemometer \iaa of the constant-temperature type 
and is described in detail in reference 15. The turbulence intensities 
were measured with an average-square computer described in reference 15. 
The correlation coefficients were measirred with a magnetic-tape correla- 
, tion computer. The anemometry eqvilpment and techniques are discussed in 

^ appendix B. 


Hertum Injection and Sampling 

Helium was supplied to the simulated point source from a standard 
14-liter bottle and a pressure-regulator system. The heliimi mass flow was 
metered rtth calibrated critical-flow orifices and injected into the 
stream through a simiilated point sovirce shown in figure 2. The injector 
was carefully streamlined to minimize flow disturbances . 

The tip of the sampling probe had a 0.040-lnch inside diameter and 
a 0.060-inch outside diameter and could be moved across the duct id.th a 
probe actuator. The probe could be positioned at any of the instrumenta- 
tion stations sho-tm in figure 1. Ohe sample was drawn from the stream 
vd.th a vacuima pump^ metered with calibrated critical-flow orifices, and 
drawn past a standard leah connected to a helium-leak detector. The leak 
detector was a msiss- spectrometer type and had a linear sensitivity to 
helium concentration. 


PROCEDURE 

Mean Stream Conditions 

Turbulence and hellvun-diffusion measurements trere made at five 
stream velocities; 50, 75, 100, 122, and 164 feet per second. The air 
tempemture varied from 70° to 78° F, and the duct static pressure varied 
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from 26.6 "to 29.0 indies of mercury absolute. These variations were 
caused by day-to-day changes in room tenrperature and barometric pressure 
and by ducting pressure drops. 


Turbulence Measurements 

The turbulence field was investigated in detail by making traverse 
measurements of the turbulence parameters at each of the Instrumentation 
stations shoim in figure 1. The central core of fully developed turbu- 
lent pipe flow should approach homogeneity and isotropy. The approach 
to these conditions wsis investigated, in addition to the turbulence 
measurements required for comparison of the Eulerian and Lagrangian cor- 
relation coefficients. The data reduction methods for the turbulence 
measurements are discussed in appendix B. 


Diffusion Measurements 

Helium- concentration profiles were measured at a series of stations 
downstream from the injector for each of the five stream velocities. In 
order to minimize turbulence effects at the Injector, most of the diffu- 
sion data were taken with the helium injection velocity equal to the 
stream velocity. With this condition, the helium concentration was too 
low for satisfactory detection at the stations farthest downstream. To 
alleviate this condition, injection flow rates were doubled for measure- 
ments at these stations. Since the relative concentration C/c^ at the 
intermediate stations was unchanged by this procedure, it was concluded 
that injection at twice stream velocity had no effect on evaluation of 
the spreading coefficient o at the stations far downstream. 

Since the concentration profiles were three dimensional, it was 
necessary to locate the peak of the profile carefully. This was done 
by rotating the sampling probe to change the vertical position of the 
probe tip and then moving the probe horizontally until the peak was found. 

The total- spreading coefficient £D was calculated from the slope 
of a strai^t line faired through each set of concentration-profile data, 

where the data were plotted an log (c/c^) against r^. The molecular- 
spreading coefficient was calculated from the relation 

= td (24) 

where the molecular-diffusion coefficient d was calculated from the 
theoretical relations given in reference 16. The turbulent- spreading 
coefficient ©51 was then calculated from equation ( 8 ) . 
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RESUIffS 

General Stream Ccaadltions 

Pitot static prole and hot-wire-anemometer measiorements showed that 
in the field of interest, the mean stream velocity had a flat profile in 
both the axial and radial directions. Ihe slight variations in the tem- 
perature (70° to 78° F) and pressure (26.6 to 29.0 in. Hg ahs) affected 
O) the molecular-diffusion coefficient of helium somewhat, and this effect 

S was included in the evaluation of the turbulent- spreading coefficient. 


Lateral Turbulence-Intensity Measurements 

The lateral turbulence intensity was measured throu^out the 

field of Interest for each of the stream velocities used. Typical radial 

and axial profiles of are shown in figures 3(a) and (b), respec- 

tively. The profiles are essentially flat in both directions. 


A variation in lateral turbiilence intensity with stream velocity 
was found and is shown in figure 4. This variation might be due to the 
turbulent field not being fully developed, since the ratio of duct length 
to diameter was only 36 in contrast to the generally accepted value of 
40 to 50. 


Longitudinal Turbulence-Intensity Measurements 

The longitudinal turbulence intensity V? was measured throu^out 
the field of interest and over the range of stream velocities. Typical 
radial and axial profiles are shOT-m in figures 5(a) and (b), respectively, 
and are essentially flat. 


Radial Measurements of the Eulerian Correlation Coefficient g 

Measurements of the Eulerian correlation coefficient g were made 
throu^out the field of interest at each of the stream velocities used. 
These measurements \rere made in the radial direction with the X-wire- 
anemometer probe, so that g was oriented as shown in the foUovTing 
sketch: 
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A typical series of g measurements along the duct axis at various 
stations do^mstream from the helium injector is shown in figure 6(a). 

The correlation-coefficient curves obtained hy fairing lines through the 
g data as shown in figure 6(a) are summarized in figure 6(h) where a 
distinct variation in the correlation coefficient with stream velocity 
can he seen. This may he due to undeveloped turbulent flow^ as is a 


similar trend of 


Vv^. 


Axial Measurements of the Eulerian Correlation Coefficient f 

The Eulerian correlation coefficient f was measured throu^out 
the field at each of the five stream velocities. These measurements 
were made in the axial direction with the single-wire-anemcmeter probe, 
so that f was oriented as follows: 


U 



I— 


A typical series of f measurements along the duct axis at various 
stations downsti’eam from the heliimi injector is shown in figure 7(a). 
The solid line in this figure is discussed in the next section. A sum- 
mary of the measured f correlation coefficients is given in figure 
7(h) and shows a variation with stream velocity similar to that of the 
g correlation coefficients in figure 6(h) . 


Turhulent-Spreadi ng-Coeff icient Measurements 

Helium- concentration profiles were measured at a series of stations 
downstream from the helium injector for each of the five stream veloci- 
ties. A typical, concentration profile is shown in figure 8, where the 
solid line represents eq.uation (12) fitted to the data hy the methods 
described in the PROCEDUEE section. As can he seen from figure 8, the 
diffusion data follow equation (12) accurately over a range sufficient 
for evaluation of the total.- spreading coefficient co. 


The turbulent-spreading coefficient was determined from the 

helium- diffusion data as described in the PROCEDUEE section. A summary 
of the turhulent-spreading-coefficient data is plotted in figure 9 as a 
function of the distance downstream from the helium injector for each 




of the five stream velocities. If the turbulence intensities 

and V^/u and the Lagrangian correlation coefficient dt were 
to remain independent of the stream velocity U, then the Taylor theory 
of diffusion predicts that a3 1 the oc^ data shown in figure 9 should 
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fall on the same curve -when plotted against the distance Xo Since both 
the turbulence-intensity and Eulerian correlation coefficients were de- 
pendent on stream velocity, the spread of the oc^ curves in figure 9 is 

to be expected. An anornalous effect is noted in figure 9 with respect 
to the (Dp curve for the stream velocity of 122 feet per second. To be 

consistent with the turbulence measurements and the other four 
curves, the cc^ in question should fall below the cDp for the stream 

velocity of 164 feet per second, but it does not. This inconsistent re- 
sult is attributed to experimental error in the diffusion data for the 
122-foot-per-second stream velocity. 


DISCUSSION 

Homogeneity of the Turbulence Field 

With the exception of two data points, the turbulent- spreading co- 
efficient was evaluated from helitnn concentrations measured within 

a 2 -inch- diameter central core of the pipe flow. This 2-inch core was 
considered to form the boundary of the turbulence field of interest. 

The consistency of the values of the turbulence parameters in this field 
is a measure of the degree of homogeneity of the field. Since the dif- 
fusion took place in the radial direction, variations in the lateral 

turbulence intensity were deemed important and are listed in the 

following table for each of the five stream velocities used: 


Stream 

velocity, 

ft/sec 

Deviation from the mean value 

of A/v in core of interest, 
percent 

50 

420.6 

75 

±8.2 

100 

410.5 

122 

47.2 

164 

414.6 


By inspection of data such as those in figure 3(a), some of the devia- 
tion shown in this table is attributed to errors in measTorement . In 

general, the radial profiles of showed sli^t trends tovrards higher 

values at the outer boundary of the core of Interest and exhibited scme- 
•vdiat higher values in the helium-injector wake for locations a few inches 
downstream from the injector. 

The correlation coefficients f and g shoired no definite trend 
In either the radial or lateral direction in the core of interest for 
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any particTilar stream velocity. The slight effect noted in the in- 
jector wake was attributed to the periodicity of the wake turbulence. 
Q?his effect decreased the value of the correlation coefficients near 
C= 0 and persisted through most of the rest of the correlation curve 
in the form of ripples superln^iosed on the generally smooth curve. 


Isotropy of the Turbxilence Field 


Isotropy is defined as the condition wherein the turbulence in- 
tensity and correlation tensor are invariant to rotation or to reflection 
in any plane throu^ the coordinate-system origin (ref. 2) . A measure 
of Isotropy is therefore found in the degree of agreement between the 
values of the turbulence intensity in the rectangular directions. Since 


only 


4 ^ 


4 ^ _and 

= 4 ^. 


obtain 
the following table: 


were measirred, radial s; 
The agreement between 



etry must be assmed to 


and 




is shown in 


ft/sec 

#/u 

yu^/u 


50 

0.0345 

0.0305 

1.13 

75 

.0340 

.0310 

1.10 

123 

.0315 

.0315 

1.00 


As discussed in the KESDITS section, the turbiilent pipe flow may not have 
been fully developed because of insufficient inlet length. This might 
explain the departure from Isotropy at the lower stream velocities shown 
in the preceding table. 

As a consequence of isotropy, a relation exists between the Eulerian 
correlation coefficients f and g (eq. (14) ) . This relation can be 
Integrated as follows to provide a numerical relation between f and- g: 

f Cg « (25) 


The agreement of the measured f and g with equation (25) was found 
by calculating an f from the measxnred g by graphical integration. 

The results of a typical calculation are shown in figure 7(a), where the 
ceilculated f is plotted for comparison with the measured f . The 
sigreement between the calculated and measured f correlation coeffi- 
cients (figs. 10 and 7(b), respectively) was best at the higher stream 

velocities, which is consistent with the agreement between and 

shown in the preceding table. 
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In sutmnftry, It appears tiiat the condition of isotropy is approached 
closely^ especially at the hi gh e-r stream velocities . The greater de- 
parture from isotropy at the lov/er stream velocities is of consequence 
in discussion later in this section. 


Agreement of the Diffusion Data with the Taylor Theory of Diffusion 

As mentioned in the INTRODUCSDION^ there are conflict in g data about 
the turbulent- diffusion process for small times (short distances down- 
stream from the diffusion source). By inspection, equation (22) indicates 
that for RTiial 1 values of t (where « l) the turbulent-spread in g coef- 
ficient CDp should be 

a>p = i v^ t^ ( 26 ) 

After correction for hot-wire length, the anemcmetry measurements of v^ 
in reference 17 were foimd to accurately predict the turbulent diffusion 
of heat from a line source by equation (26). In the similar experiments 
reported in references 7 and 8, the turbulent- spreading coefficients of 
the heat \rake -iTere greater than those predicted from anemometry measure- 
ments and equation (26) . 

The pulse-like natixre of the instantaneous temperature measured 
close be hin d the line source of heat indicated that the heat wake was a 
randomly raving sheet having a finite thickness due to molecular diffu- 
sion (refs. 7 and 8). Reference 8 suggests that the smallest eddies of 
the tiarbulence distort the sheet by rotation and strain so that the 
spreading coefficient o:^ predicted from the anemometer measurements by 

equation (26) is less than the cnp found by temperature-distribution 

measurements. The criterion for this effect is that only eddies smaller 
than the molecular wake will contribute to the distortion. 

The criterion of microscale sTnall ccnrpared .with molecTolar-wake width 
treis evaluated in the present investi^tion by determining the microscale 
from the correlation coefficient f and calculating the molecular- 

wake width as described in the PROCEDURE section. The microscale 
was determined from the definition given in reference 18 



11m ^ - f\ 

S-oUs ) 


(27) 


and was found to be approx±mately 0.015 foot for all stream conditions 
used in the investigation. The maxlnTum molecular spreading ras en- 
countered at t = 0.0497 second, where the standard deviation of the 
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molecular wake was calculated to te 0.0088 foot. These calcula- 

tions indicate that the effect described in reference 8 did not occiar 
in the present investigation, since the mayimimi molecular- wake width 
was less than the microscale of the turbulence. This conclusion is 
borne out by a ccmparison between the determined from the heli\mi- 

diffusion data and those calculated from the anemometry meastirements 
and equation (26) . This coDg)arison, which shows acceptable agreement 
between the anemometry data and those measured by diffTision, is tabu- 
lated as folJLows: 


Stream 

velocity, 

ft/sec 

Time, 

■b, 

sec 

Spreading coefficient, 
sq ft 


(from 

diffusion 

data) 

(calculated 

from 

anemometry 

data) 

0^ 

50 

0.003 

0.00001263 

0.0000134 

1.06 

75 

.002 

.00001134 

.000013 

1.15 

122 

.00123 

.00001153 

.0000111 

.96 

164 

.000937 

.00001063 

.0000106 

1.00 


Comparison of the Shapes of the Lagranglan and Eulerlan 
Correlation Coefficients 

As dlscxissed in the THEORETICAL AMLISIS section, a ccraparlson of 
the shapes of the Lagranglan and Eulerlan correlation coefficients may 
be made by comparing the spreading coefficient odqi with an 0 )p calcu- 
lated from the Eulerlan .correlation coefficient by equation ( 23) . The 
Oij. calculated from the f correlation coefficients shown in figure 10 

eire shown in figure 11 as a function of the separation distance K. 

The turbvilent-spreading coefficients oorii measured by helium diffusion 
and shorn in figure 9, are plotted against time t in figure 12 by us- 
ing the transformation x = Ut. In order to determine the validity of 
equation (2l), the separation distance C from figure 11 and the time 
t from figure 12 were cross-plotted (fig. 13) for equal values of o:^ 

and ccqi. Equation (2l) liiq)lies a linear relation between the values 
of € and t corresponding to equal values of oqj. and o^. The ex- 
perimental points on figure 13 can be approximated by strai^t lines, 
as shown, and equation (2l) is thiis substantiated. 
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The factor B in equation (21) was evaluated from the slopes of 
the faired lines in figure 13 by using values of the lateral turbulence 

intensity from figure 4. The factor B is plotted in figiire 14 

for each of the five stream velocities. 


Additional Valuation of the Proportionality Factor B 

As a part of another research program, concentration profiles of 
acetone vapor were measured downstream from a simulated point source in 
a 6 - inch- diameter duct. Turbulent- spreading coefficients were deter- 
mined from the acetone-vapor concentration profiles at distances of 20.3 
and 24.8 inches downstream from the point of Injection. 

Longitudinal turbulence intensities were measured throu^out 

the duct and had approximately flat radial and axial profiles in the 

field of interest. The turbulence intensity aP had a value of 14 
feet per second at the stream conditions listed in figure 15. 

The correlation coefficient f -^fas also measured at the stream 
conditions listed above, along the duct axis. These measurements are 
plotted in figure 15 for a series of downstream stations and shovr little 
c h an g e with position in the tixrbulence field. 

On the basis of measurements of longitudinal intensity and the f 
correlation coefficient, values of a>f were calculated from equation 

(23) for comparison with the values obtained from the acetone-vapor 

diffusion experiment. The cross plot that was made between the two sets 
of 0 ^ and cdj values gave two values of C for respective values of 

t. The (^,t) coordinates fell on a stral^t line having a slope such 
that the factor B in equation (21) had a value of 0.62. This value 
of B agrees well with those shown in figure 14. 

In summary, the data show that the Lagrangian and Eulerlan correla- 
tion coefficients have similar shapes connected by an eq,uatlon such as 
equation (21) and that the factor of proportionality in this equation is 
roughly 0.6. 


CONCLUSIONS 

1. The Lagrangian and Eulerian correlation coefficients in homo- 
geneous, isotropic turbulence have similar shapes. 
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2. The shapes of the Lagraogian and Eulerlan correlation coeffi- 
cients i-rere found to he related hy a linear proportionality between 
their time and space coordinates as shown hy the expression 

f(?) =3if — i=\. The proportionality factor B is roughly 0.6 throu^- 
out a range of tirrhulence Intensities from 1.8 to I 4 feet per second. 

3. The linear relation hetvreen the shapes of the Lagrangian and 
Eulerian correlation coefficients provides a means for the quantitative 
solution of practical mixing problems from hot-wire-anemometiy data. 


Levri-s FI 1 ght Propulsion Laboratory 

National Advlsoiy Committee for Aeronautics 
Cleveland^ Ohio, August 5, 1955 
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APPMDIX A 
SICMBOLS 

Tlie foll04d.ng symbols are used in tMs report: 

B factor of proportionality in linear relation betireen shapes of 

Lagrangian and Eulerlan correlation coefficients 

C concentration of diffusing gas or vapor 

D turbulent diffusion coefficient, sq. ft/sec 

d molecular diffusion coefficient, sq. ft/sec; or cylinder diam- 

eter, ft 

E direct current voltage, v 

root-mean-sq.uare of alternating voltage, v 

f,g Eulerlan double velocity correlation coefficients defined in 
manner of von KarmfCn and Hoi-rarth 

G electronic amplification factor 

K constant correcting for division of current in .Wheatstone bridge 

L Eulerlan scale of turbulence 

-Sf’ Lagrangian scale of turbulence 

n frequency of vortex- shedding 

^ Lagrangian correlation coefficient 

r radieil distance from axis in cylindrical coordinates, ft 

S Strouhal number 

t time, sec 

ij|c 

t characteristic time of turbulence, sec 

U stream velocity, ft/sec 

u,v,w turbulent-velocity vectors in x, y, and z directions, respec- 
tively, ft/sec 
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W instaataneoiis point- soiorce strength, lb 

■w injection rate, Ib/sec 

X distance in axial direction in cylindrical coordinates, in. or 

ft 

Y fluid particle displacement, ft 

a angular displacement of X-wire ariay -with respect to stream 

direction, radians 

X microscale of turhulence, ft 

K separation distance "between velocity vectors in Eulerian douhle 

velocity correlation coefficients , ft 

2 electrical resistance, ohms 

CD spreading coefficient, sg. ft 

Subscripts : 

h voltage across Wheatstone bridge 

c cylinder diameter 

D difference between voltages across two Wheatstone bridges 

d molecular 

f pertaining to correlation coefficient f 

i summation convention 

o peak, or no-flcjw condition 

T turbulent 

t time 

g at distance g 

Baonred symbols Indicate time average 
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APPENDIX B 


HOT-WIRE ANEMCMSEEX 
Anemametry Equipment 

THe hot-wire-anemometry equipment is shown diagraimnatically in fig- 
ure 16. The aiiq)lifiers were of the constant-tenqierature type with a 
frequency-response flat within 1 decihel from 0 to 100,000 cycles per 
second. The bridge and anq)lifler design is discussed in detail in ref- 
erence 15. 

The turbulence intensities were measured with the average-square 
ccmq)uter described in reference 15. The computer had a frequency- 
response -curve flat within 0.3 decibel from 20 to 10,000 cycles per 
second, a range of 0.02 to 5 volts rms, and a variable- capacitance 
meter damping circuit. The conq)uter was calibrated with a sine wave, 
but measured true root-mean- square voltage of complex waves. 

T?he difference amplifier was a standard design with a gain of 1.0 
and a flat frequency- response curve within 0.5 decibel between 20 and 
20,000 cycles per second. Its function was to obtain the Instantaneous 
difference between the complex voltages from two hot wires. 

The magnetic-tape recorder was of the hi^-speed two-channel type 
and was used to store the complex voltage signals from the hot-wire 
anemometers for later analysis. 

The correlation computer was used to determine the correlation co- 
efficients from the anemometer voltage signals stored on the m ag n etic 
tape. The curve for over-Eill frequency response of the computer was 
flat within 0.5 decibel from 25 to 5000 cycles per second. The operation 
of the computer is described later in this appendix. 

The X-^/lre- and single-wlre-anemcmeter probes for measuring the 
lateral and longitudinal turbulent velocities, respectively, are shorn 
in figure 17. 


Turbulence-Intensity Measurement 

The longitudinal turbulence intensity was determined from the equa- 
tion 
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where and Ejj ^ are direct- current voltages measured at flow 

and no-flow conditions, respectively, and *1^ is the root-mean-square 
voltage given hy the average-squane computer. Equation (B1) Is based 
on King's equation and assumes that the natural convection and com- 
pressibility effects on the forced convection frcan the hot wire are 
negligible . 


The lateral ttirbulence intensity ^"^2 determined by the method 

developed in reference 19 and discussed in reference 20 that relates the 
lateral velocity fluctuation to the voltage difference between the two 
hot wires in the X-array through the equation 



Where the root-mean-square voltage >^E^ is an instantaneous difference 
voltage read from the average- squaxe coiDputer, and the term (dE^)/(da) 
is the slope of the direct-current difference voltage against the angu- 
lar position of the X-array with respect to the stream- velocity vector. 


Correlation- Coefficient Measurement 

The correlation coefficients were measured by a method introduced 
in reference 21 that inserts an artificial time lag in the tTorbulence 
signal from a single anemometer probe. This method has definite op- 
erational advantages over the conventional two-probe method, in that 
only one stationary probe is required and that the anemometer signal 
may be stored for analysis at convenience. 

The correlation computer used in the present investigation, which 
was designed by D. F. Berg and C. C. Conger of the. Lewis laboratory, is 
shown dlagrammatically in figure 18. The magnetic tape supplied to the 
computer had identical signals on both channels. The computer \ms 
equipped with one fixed and one moveable play-back head, so that an ad- 
justable time delay could be Introduced between the otherwise identical 
signals passing through the two channels of the computer. The out-of- 
phase signals were amplified in their respective channels by anpliflers 
matched to the high-speed tape recorder. After amplification, the two 
signals were multiplied instantaneously by an electronic analog multi- 
plier. The analog multiplier was of the modulated carrier-vrave type 
with a special wave foim having a freq.uency of 2 megacycles. The time 
average value of the product of the two signals was a correlation 
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coefficient which related the velocities of two different particles at 
two different times. The time intenral between the two velocities was 
obtained by dividing the separation distance between the two pickup 
heads by the tape speed. The space interval between the two velocities 
was obtained by multiplying the time interval by the stream velocity. 

Since the turbulence intensity was much less than nnity, 

the lateral displacement of fluid particles during the time lag Interval 
t of the con 5 )uter was smal 1 compared with the longitudinal, space in- 
terval given by 5 = Ut. From this consideration, the correlation co- 
efficient given by the ccmputer should approximate the true Eulerlan 
coefficient. This concl\ision is substantiated by measurements reported 
in reference 22 ^dilch show good agreement between correlation coeffi- 
cients obtained from the computer and a single probe and those obtained 
by the conventional two-probe method. 
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Section A-A 


Figure 2. - Helium injector. 
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Distance dcjvmstreeun from helium injectoo", x.j in. 

(b) Axial profile along duct axis. Stream velocity, 75 feet per second. 
Figure 5. - Typical longitudinal turbulence intensity profiles. 
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Figure 7. - Axial meetaurementB of Eulerlan oorrelation coefficient f. 





Radial difftance from conceirtratLon profile center, r, in. 

Plgure 8. - Typical radial profile of helium concentration. Stream velocity, 50 feet per aecond; dlatanca 
dovnatream from injector, 12 Incbee. 
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Distance downstream from point source, x, ft 


Figure 9 . - Summary of turDulent-spreading-coefficient data ob- 
tained from helium concentration profiles . 
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rjLgura 10. - Sumnary of oorrelatlon coefficient f calculated from neasm^ed correlation coefficient 
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Figure 12. - 3mninary of tine variation of turbulent-Bpreadlng ooofflolant ineaBured by hallum dlffuBlon. 
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(b) Circuit used for longitudinal turbulence measurements 
Figure 16. - Hot-wire-anemometer circuits. 
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Tungsten array 

Dlam., 0.0002", 
length, 0.080" 



j. 


Air f Ixjw 



dlam. tapered 
from 0.031" to 
0.007" at tips 


(a) X-wlre-anemometer prohe. 


wire 



■0.63"* 





20 " 


length, 0.080' - 1 aigjports, 

I diam. tapered 
from 0.031" to 
0.007" at tips 


‘-Diam., 

0.125" 


^Diam., 0.25" 


2 


(h) Single-wire-anemometer prohe. 
Figure 17 . - Hot-wire-anemometer prohes 
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Figure 18. - Block diagram of correLation computer 
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